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Exact solution of the linearized equations for steady-state transport in semiconductors yields two 
modes that vary exponentially in space, one involving screening (without entropy production) and 
one involving diffusion and recombination (with entropy production). Neither mode is quasineu- 
tral. For constant surface photoexcitation with generation of electrons and holes, the steady-state 
response is a linear combination of these modes, subject to global electroneutrality. The resultant 
charge separation produces a voltage difference across the sample (the Dember effect). 
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The Dember effect is the voltage difference that devel- 
ops across the bulk of a material that steadily absorbs 
light at one surface. M It has been observed in insula- 
tors, ||,|| high resistivity semiconductors, |Qg and Cqq 
thin films, [pi An extensive discussion of photovoltages, 
both at surfaces and in the bulk, is given in Ref. M. 

An intuitive physical picture of the Dember effect is as 
follows. H Incident light produces equal numbers of elec- 
trons and holes at the surface. The higher mobility elec- 
trons travel further than the lower mobility holes. The 
resultant charge separation amounts to a dipole layer. 
As a consequence the illuminated surface is at the higher 
voltage, in agreement with experiment. However, con- 
ventional theories of the Dember effect, [[7||| which as- 
sume quasineutrality (local electroneutrality), yield no 
dipole layer; moreover, they yield a non-zero net charge. 

The present work obtains exact solutions for the sur- 
face modes of the linearized macroscopic transport equa- 
tions. It differs from previous analytic work in that: (1) 
It explicitly satisfies the principles of irreversible thermo- 
dynamics, which constrain the form for the charge-carrier 
fluxes and recombination rates by the condition that the 
rate of entropy production be non-negative; (2) In solv- 
ing the transport equations, it employs both electron and 
hole densities n and p as variables, rather than assum- 
ing deviations from equilibrium satisfy quasineutrality 
(Sn w Sp); (3) Using the resultant two steady-state sur- 
face modes, it explicitly enforces overall electroneutrality. 

One surface mode corresponds to ordinary screening by 
a multi-charge-carrier system (i.e. Debye-Hiickel screen- 
ing), without fluxes or entropy production: the screening 
mode. The other surface mode has both electron and 
hole flux as well as entropy production: the diffusion- 
recombination mode. The characteristic length I of the 
screening mode typically is shorter than the characteris- 
tic length L for the diffusion-recombination mode. Nei- 
ther mode is locally electroneutral. Applied to the Dem- 
ber effect, the light generates a diffusion-recombination 
mode with a relatively extended negative charge and a 
screening mode with a relatively compact positive charge, 
in agreement with the physical picture given above. 

Irreversible thermodynamics. Consider a uniform 
semiconductor with an ideal surface having no extrinsic 



surface states and no charged intrinsic states. The re- 
combination rate r is the same for both electrons and 
holes. With u the energy density, T the temperature, s 
the entropy density, and fi e and jlh the electron and hole 
electrochemical potentials, the fundamental thermody- 
namic differential for this system is 



du = Tds + fi e dn + jlhdp. 



(1) 



With n e and \ih the chemical potentials (to be distin- 
guished from the diffusivities /j, n and fx p ), 

p, e =fi e ~ e<j), jlh = Hh + ecj). (2) 

Here the electrical potential cf> satisfies Poisson's equation 

V 2 0=-^ = - — (P-U + N+-N-), (3) 

where e is the semiconductor dielectric constant, eo is 
the permittivity of free space, and the charge density 
p — e(p — n + N^[ — N~). N^ and N~ are the respective 
densities of ionized donors and acceptors, and jl e and jlh 
are often called quasi-Fermi energies. 

The conservation laws for this system are 

9 t u+a i j? = o, d tS + dij! = ^>o, (4) 

ftn + fitf? = r, d tP + dijf = r. (5) 

Here jf is the energy flux density, jf is the entropy flux 
density, V is the local rate of heat production (V/T is 
local rate of entropy production), j" is the electron num- 
ber flux density, and jf is the hole number flux density. 
The fluxes, r, and V are to be determined. 

The time-derivatives of (Q) and (@)-(|§) lead to 

< V = -diJT + TdiJt -Mr- diJT) - Mr - djf ) 

= -di(jr-TJt-M?-hJ!) 

-JidiT - j?d^ e - jfd^h - r(fl e + iX h ). (6) 

Expressing V as a non-negative quadratic form re- 
quires that jf = Tjf + jl e jf + jihjf; that 
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where n > is the thermal conductivity; that 



xp _ 






(8) 



where a nn > 0, a pp > 0, and a nn a pp > a np a pn (the a's 
are related to the diffusivities, and a np = a pn by the 
Onsager symmetry principle |9j]); and that 

r = -\(£i e + [i h ) = -\(fi e + pt h ). (9) 

where A > is related to the electron and hole lifetimes 
T n and r p . In equilibrium r = 0, so (|9j) implies the well- 
known results that fi e = —fih and \i e = — /i/, in equilib- 
rium. We now express r, jf, and jf in more conventional 
form. 

Linearizing (|9j) about equilibrium with Su e = 
(du e /dn)5n and Sn^ = (d "fih / 'dp)5p yields 



dn dp 



X 



ftT' T " = V- (10) 



For a dilute semiconductor, as considered here, 
(dn e /dn) ss k B T/n and (du h /dp) « k B T/p , so r„ = 
(n Q /k B TX) and r p = (p /k B TX). 

Now set a„ p = a p „ w 0, a„„ = D n (dn/du e ), and 
a pp = D p {dp/d[ih), where £>„ and D p are the electron 
and hole diffusivities. Next, define the electron and hole 
mobilities /x„ and u pi which satisfy the Einstein relations: 



eD n dn eD 7t 



Hn = 



_eDp&p_ eD p 
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n <9/x e fcflT' ' " p dfih k B T 
Then in one-dimension jf and jf become pJLfil 



J., 



-D n d x n + n n nd x (j), 



Jl = -D p d x p - fi p pd x cj). 
Note that Poisson's equation linearizes to 



(12) 
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(13) 



Static surface solutions in one-dimension. One 

steady-state solution that automatically satisfies (J5J) , (g) , 
and (||) has <5/2 e = Sfih = and jf — jf — 0. Thus the 
system is in local equilibrium, with no entropy produc- 
tion. This corresponds to ordinary screening, and thus 
we call it the screening mode, with subscript S. From 
(O), its potential 54>s satisfies 

V 2 ^ s = —{£- + ^)Hs = q 2 s 6<f> s , (14) 

e e dp h d/i e 

where 
22 ,2 2 _ e2 9n 2 _ e 2 dp ,-. K s 

IS — ISn + IShi QSn = ~~^~~' ^Sp = ~~~JT~ • ( 15 J 

eoe clfi e r e e Ofi h 



For the present system, q 2 Sn — (e 2 n /e ek B T), q 2 s = 
(e 2 p /e ek B T), so q 2 s = (e 2 /e ek B T)(n +p ). The solu- 



(16) 



tion to (|14|) that goes to zero as x — > cxj is 
50s = A s exp(-g s a;). 



I = q s is the screening length. From (16) and dlq), the 
screening mode has charge density 

Sps = -eoeV 2 0s = -{e Q e)q 2 s A s exp (-qsx). (17) 

Here <5ns = (dn/du e )Su e = {eno/k B T)5(j)s and (5ps = 
(dn/d/j, h )5u h = -(en /k B T)S(f>s, so 6n s /6p s = -n /p . 
The second steady-state solution to (||), (g), (g) and 
(O) has (S/i e and 5fih non-zero, so the fluxes jf and jf 
are non-zero. With a np — a pn — 0, combining (ph, (H) 
and (g) yields 



- a nn dl8jl e = -X(8jl e + Sflh), 
-a pp d 2 x 5jl h = -X(Sp, e + Sfih). (18) 

Let us denote properties of this diffusion-recombination 
mode by the subscript D. Now set 
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Then, with 

5 fi e = ^4_Dn exp (-q B x), Sflh = A Bp exp {-q B x), (20) 
© yields 
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and A Dp /q 2 Dp = A Dn /q 2 Dn . 

The corresponding electrostatic potential 
the form 



to = A r ,exp(-g £) a;), 



where 
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(21) 
i>D takes 

(22) 

(23) 



L = qfj is the diffusion-recombination length, normally 
called the diffusion length. By (22) and (3), the diffusion- 
recombination mode possesses charge density 



SpD = -(e e)q 2 D A D exp {-q D x) 



(24) 



Moreover, Sn D /Sp D = (n /p )(q 2 Sp - toJ/(<?L - top)- 
This is independent of no for fixed noPo = n f ■ 

Surface mode amplitudes. Let light incident at the 
surface x = produce equal electron and hole fluxes G, so 
jf = jf = G at x = 0. On setting a np — a pn — 0, (|J) and 
(|20j) yield the diffusion-recombination mode amplitudes 
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(25) 



Overall electroneutrality, or 
5pd(x)), implies that 

qs 



f™dx(5p s (x) 



From (17) and (C 
density is 



(26) 



and from (pS-Ea) , the total charge 
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where the mobilities have been employed. Typically p n > 
Up, so p(x — 0) is positive, as expected if the higher 
mobility charge-carrier preferentially leaves the vicinity 
of the surface. Note that Sp(x) changes in sign as x 
increases, as needed to produce a dipole layer. 

From (|l7|) and (|24|), the total electrostatic potential is 



6cf>(x) 



Ge 2 (p n - pp) ggV 
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(28) 



To compare with experiment, (pq) gives the Dember 
voltage 4>Dem = d<f)(x = 0) as 



-'Dem 



Ge 2 (p n - pp) 



1 



e ek B T p n p p q s qD qs + to 



(29) 



Discussion. Eq.(29) makes numerous predictions. 
First, since typically p n > p p , <pDem is positive, as ex- 
pected for a dipole layer with positive charge closer to 
the surface. Further, since typically qu <C qs, the term 
9sto('Zs' + to) ~ ?Ito- The dependences of qs on n and 
e, and of qn on p and r, show that <pDem varies inversely 
with carrier density n, as the square root of the charac- 
teristic recombination time r, as the inverse square root 
of the characteristic mobility p, and is independent of e. 

The dependence on n explains why <pDem is observed 
in materials with relatively low carrier densities. The 
dependence on r indicates that 4>Dem is larger for longer 
recombination times; the dependence on p indicates that 
4>Dem is larger for lower mobilities. 

Figures 1-4 are computed for an intrinsic semiconduc- 
tor with equilibrium carrier number density no — Po — 



10 16 m 3 , recombination time r 



1(T 



elec- 



tron mobility p n = 1000 cm 2 /(V-s), hole mobility p p = 
200 cm 2 /(V-s), room temperature k B T — 0.0253 eV, and 
dielectric constant e = 10. These values give I — 26 A pvn 
and L = 64.8 /mi. From (28), the ratio of the contri- 
butions of the diffusion mode and the screening mode 
to the Dember voltages is q^ 1 /q^ 1 = L/l. Taking 
G = 10 17 m^s" 1 gives cj> Dem = 9.22 mV, of which 
the diffusion mode is responsible for 15.56 mV and the 
screening mode —6.34 mV. The theory, which is valid for 
5p e and 5p,h small relative to k B T, becomes invalid for 
significantly larger G. 
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FIG. 1. Photo-induced electrical potential 
and charge density 5p from (27). 
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FIG. 2. Electron and hole number densities n and p corre- 
sponding to Figure 1. 



Charge Densities of Screening and Diffusion Mode 



Charge Density for Different Mobilities 
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FIG. 3. Charge densities for the screening and diffusion 
modes, from (17) and (24). 
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FIG. 4. Charge density Sp from (27) for different hole mo- 
bilities. 



Fig.l presents Scf> and Sp. Note that S(f> is monotonic, 
whereas Sp changes sign, so that the system can be overall 
neutral. Fig. 2 presents the number densities n and p 
associated with 5p of Fig.l. The curves sketched (but 
not computed) in Ref. g (Fig.V.3) agree qualitatively 
with Figs.l and 2. 

The individual charge densities p due to the diffusion- 
recombination and screening modes are presented in 
Fig. 3. Note that Sns/Sps = —1 and Shd/Spd = 1-306 
for the present parameters. (Increasing both mobilities 
by a factor of ten would give Sn^/Spo = 1.022.) Hence 
quasineutrality is not a good approximation for either 
mode. The charge density for the screening mode is 
positive and concentrated near the surface, whereas the 
charge density for the diffusion-recombination mode is 
negative and more extended. Note that numerical solu- 
tions of the macroscopic transport equations for semicon- 
ductors normally do not assume quasineutrality, so they 
retain the full physics. Jl^] However, they do not permit 
extensive physical interpretation. 

Fig. 4 presents the bulk charge density Sp for pt n and pt p 
as above (p, p /p n — 0.2), and for three additional values 
of p p /p, n . 5p and 4>Dem increase with the difference in 
the electron and hole mobilities. 

Surface solutions also are relevant to current flow 
across interfaces. Already in the nineteenth century, 
Maxwell noted that for current flow between two single- 
carrier materials charge builds up at the interface; how- 
ever, without a more microscopic theory it is not possi- 
ble to determine how much charge is associated with each 
material. A similar effect must occur at a semiconductor- 
metal interface, even without the complexities associated 
with modification of the Schottky barrier. Surface solu- 
tions likely are present in recent numerical studies show- 
ing that spin injection in a semiconductor (e.g., by pref- 
erential absorption of polarized light) persists across a 



p — n junction. |T|] Such a system, with up and down 
electrons and with unpolarized holes, should have three 
surface modes: one screening mode and two diffusion- 
recombination modes. The rapid falloff of spin polariza- 
tion at one end of the system is evidence for surface so- 
lutions induced by what the authors of Ref. |13[ consider 
to be an artificial boundary condition. 

The present multi-carrier system, with a finite 
diffusion-recombination length L, is non- neutral only 
over a finite distance from the disturbance. On the other 
hand, a multi-carrier system with no recombination (and 
thus L — y oo ) can be non-neutral in the bulk, far from a 
disturbance. [H 
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